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METHOD FOR DETERMINING MOLECULAR PROPERTIES OF 
HYDROCARBON MIXTURES FROM NMR DATA 

Background of Invention 

Field of the Invention 

[0001] The present invention is related to the field of data processing methods for oil 

well logging and sampling. More specifically, the present invention relates to methods 
for determining properties of hydrocarbon mixtures and crude oils including molecular 
composition, molecular size, molecular weight, and molecular carbon number using 
nuclear magnetic resonance (NMR) data. 

Background Art 

[0002] Oil well logging and sampling tools include nuclear magnetic resonance (NMR) 

instruments. NMR instruments can be used to determine properties of earth formations, 
such as the fractional volume of pore space, the fractional volume of mobile fluid filling 
the pore space, and the porosity of earth formations. In addition, NMR data may be used 
to assess the content of brine and hydrocarbons in the formation. General background of 
NMR well logging is described in U.S. Patent No. 6,140,817, assigned to the assignee 
hereof. 

[0003] The signals measured by nuclear magnetic resonance (NMR) logging tools 

typically arise from the selected nuclei present in the probed volume. Because hydrogen 
nuclei are the most abundant and easily detectable, most NMR logging tools are tuned to 
detect hydrogen resonance signals (from either water or hydrocarbons). These hydrogen 
nuclei have different dynamic properties (e.g., diffusion rate and tumbling/rotation rate) 
that are dependent on their environments, such as the chemical structure and size of the 
molecules in which they reside. The different dynamic properties of these nuclei 
manifest themselves in different nuclear spin relaxation times (i.e., spin-lattice relaxation 
time (Tl) and spin-spin relaxation time (T2); spin-lattice relaxation is also referred to as 
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longitudinal relaxation, and spin-spin relaxation as transverse relaxation). For example, 
molecules in viscous oils cannot diffuse or tumble as fast as those in light oils. As a 
result, they have relatively short relaxation times. These observations suggest that NMR 
data (e.g., relaxation times) can provide information on molecular properties of 
hydrocarbons in the earth formations. 



Summary of Invention 

[0004] One aspect of the invention relates to methods for estimating molecular properties 

such as composition, size, carbon number and weight in a mixture (e.g., crude oils) from 
NMR data. A method for determining molecular properties in a mixture of hydrocarbons 
includes measuring NMR data of the mixture using an NMR tool or a laboratory NMR 
flj instrument; deriving at least one parameter for each observed constituent in the mixture 

gi from the NMR data; and calculating a molecular property for each observed constituent 

2f in the mixture from the at least one parameter. Methods according to some embodiments 

yy 

s of the invention use correlations between relaxation times and molecular properties 

u 

ry and/or between diffusion rates and molecular properties. 

U 

03 [0005] Other aspects of the invention would become apparent from the following 

B 

P description, the drawings, and the claims. 

Brief Description of Drawings 

[0006] FIG. 1 is a diagram of a nuclear magnetic resonance tool in a borehole penetrating 

earth formations. 

[0007] FIG. 2 is a diagram illustrating components of a nuclear magnetic resonance tool. 

[0008] FIG. 3 illustrates a schematic diagram of a circuitry of an NMR tool for producing 

the RF pulses and for receiving and storing the spin echoes. 

[0009] FIG. 4 is a diagram illustrating a common pulse sequence for measuring 

transverse relaxation times of NMR signals and the resultant spin echoes that can be used 
to derive transverse relaxation times of nuclear magnetic signals. 
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[0010] FIG. 5 shows a flow chart of steps involved in the methods according to 

embodiments of the invention. 

[0011] FIG. 6 shows results of estimating molecular sizes in mixtures using methods of 

the invention as compared with those from gas phase chromatography. 

Detailed Description 

[0012] FIG. 1 shows a nuclear magnetic resonance (NMR) logging tool 30 for 

investigating earth formations 31 traversed by a borehole 32. The NMR logging device 
30 is suspended in the borehole 32 on an armored cable 33, the length of which 
substantially determines the relative depth of the device 30. The cable length is controlled 
by suitable means at the surface such as a drum and winch mechanism (not shown). 
Surface equipment 7 can be of conventional type and can include a processor subsystem 
which communicates with downhole equipment including NMR logging device 30. 

[0013] The NMR logging device 30 can be any suitable nuclear magnetic resonance 

logging device; it may be one for use in wireline logging applications as shown in FIG. 1, 
or one that can be used in logging while drilling (LWD) applications. The NMR logging 
device 30 typically includes a means for producing a static magnetic field in the 
formations, and a radio frequency (RF) antenna means for producing pulses of magnetic 
field in the formations and for receiving the spin echoes from the formations. The means 
for producing a static magnetic field may comprise a permanent magnet or magnet array, 
and the RF antenna means for producing pulses of magnetic field and receiving spin 
echoes from the formations may comprise one or more RF antennas. 

[0014] A schematic representation of some of the components of an NMR logging device 

30 is illustrated in FIG. 2, which shows a first centralized magnet or magnet array 36 and 
an RF antenna 37, which may be a suitably oriented coil or coils. FIG. 2 also illustrates a 
general representation of closely-spaced cylindrical thin shells, 38-1, 38-2. . .38-N, that 
can be frequency selected in a multi-frequency logging operation. One such device is 
disclosed in U.S. Patent No. 4,710,713. In FIG. 2, another magnet or magnet array 39 is 
shown. Magnet array 39 may be used to pre-polarize the earth formation ahead of the 
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investigation region as the logging device 30 is raised in the borehole in the direction of 
arrow Z. Examples of such devices are disclosed in U.S. Patent Nos. 5,055,788 and 
3,597,681. 

[0015] FIG. 3 illustrates a schematic of a circuitry of an NMR tool for producing the RF 

pulses and for receiving and storing the spin echoes. One skilled in the art would 
appreciate that any other suitable circuitry could be used without departing from the 
scope of the invention. 

[0016] In FIG. 3, a downhole processor 210 has associated memory, timing, interfaces, 

and peripherals (not separately shown), as known in the art. The processor subsystem 210 
is coupled with telemetry circuitry 205, for communication with a processor on the 
surface (not shown). The pulse forming circuitry includes a variable frequency oscillator 
220 which, under control of processor 210, produces radio frequency (RF) signals at the 
desired frequencies. The output of oscillator 220 is coupled to a phase shifter 222 and 
then to a modulator 230, both of which are under the control of processor subsystem 210. 
The phase shifter 222 and modulator 230 can be controlled, in a manner known in the art, 
to produce the desired pulses of RF field, for example the 90 degree and 180 degree 
pulses for Carr-Purcell-Meiboom-Gill (CPMG) types of sequences or any other desired 
NMR pulse sequences. The output of modulator 230 is coupled, via a power amplifier 
235, to the RF antenna 240. A Q-switch 250 can be provided to damp the RF antenna 
system to reduce antenna ringing. The antenna 240 is also coupled with a receiver 
section via duplexer 265, the output of which is coupled to receiver amplifier 270. The 
duplexer 265 protects the receiver amplifier 270 from the high power pulses which pass 
to the RF antenna 240 during the transmitting and damping modes. During the receiving 
mode, the duplexer 265 acts as a low impedance connection from antenna 240 to the 
receiver amplifier 270. The output of receiver amplifier 270 is coupled to a dual phase- 
sensitive detector 275, which also receives, as a reference, a signal derived from the 
oscillator signal. The detected output is coupled to analog-to-digital converter 280, the 
output of which is a digital version of the received nuclear magnetic resonance signal. 
Although the logging device or tool 30 is shown as a single body in FIG. 1, it may 
alternatively comprise separate components, and the tool may be combinable with other 
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logging tools. Also, while a wireline is illustrated, alternative forms of physical support 
and communicating link can be used, for example in a measurement while drilling 
system. 

[0017] Several NMR parameters may be measured that can be used to derive formation 

properties. Most NMR logging operations measure the spin-lattice (longitudinal) 
relaxation times (Tl) and/or spin-spin (transverse) relaxation times (T2) of hydrogen 
nuclei. In addition, some NMR logging tools may provide a ratio of T1/T2 directly, and 
other NMR tools may provide diffusion constants (D). These NMR data (Tl, T2, T1/T2, 
and D) are all applicable to the embodiments of the present invention, though the 
M; following discussion uses T2 relaxation times to illustrate the present invention. 

s 

|j [0018] Various pulse sequences are available for measuring the NMR relaxation times, 

dy For example, Tl relaxation may be measured using an inversion-recovery or a simple 

U spin-echo pulse sequence or any of their derivatives. The T2 relaxation is often measured 

M from a train of spin-echoes that are generated with a series of pulses such as the Carr- 

s 

U Purcell-Meiboom-Gill (CPMG) pulse sequence or some variant of this. The CPMG pulse 

p; sequence is well known in the art. (See Meiboom, S., Gill, D., 1958, "Modified Spin 

M Echo Method for Measuring Nuclear Relaxation Times," Review of Scientific 

Q 

U Instruments, 29, 688-91). As illustrated in FIG. 4, the CPMG pulse sequence generates a 

train of spin echoes, whose amplitudes exponentially decay as a function of time. The 
exponential decay life time is referred to as a transverse relaxation time, T2. Thus, T2 
measurements are accomplished by analyzing the amplitudes of spin echoes thus 
obtained. 

[0019] As shown in FIG. 4, in a CPMG sequence, the first RF pulse applied to antenna 

(37 in FIG. 2) is a 90-degree pulse, which reorients the hydrogen nuclei onto a plane 
("transverse plane") perpendicular to the static magnetic field produced by the permanent 
magnet 36. Shortly after the initial 90-degree pulse, a train of 180-degree pulses (with a 
delay time between the successive 180-degree pulses, Ti 8 o, approximately twice as long 
as the initial delay between the 90-degree and the first 180-degree pulses, T 90 ) is applied 
to the antenna (37 in FIG. 2). Each of these 180-degree pulses results in a spin echo - a 
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growth and subsequent decay of the detected signal magnitudes. During these 
measurements, the nuclear spins in the transverse plane gradually decrease amplitudes 
due to spin-spin interaction and other relaxation mechanisms. Consequently, each 
successive spin-echo has a lower amplitude than that of the preceding one. T2 relaxation 
time (the transverse relaxation time) information is then derived from analysis of the 
exponential decay profile. 

[0020] Once NMR data (e.g., Tl, T2 relaxation times, T1/T2 ratio, or diffusion rates) are 

collected, they are analyzed with an inversion method to derive the earth formation 
information. Any of the inversion methods known in the art are suitable. For example, 
U.S. Patent No. 5,291,137, issued to Freedman and assigned to the same assignee hereof, 
discloses a "windows" processing method. This "window" processing method is suitable 
for most NMR data analysis. 

[0021] Transverse (T2) relaxation in liquid is mainly through dipole-dipole interactions, 

which are influenced by the dynamic properties of the molecules (e.g., diffusion rates and 
molecular tumbling rates) and the fluids (e.g., viscosity). Thus, NMR data (especially, 
T2) may be used to provide information on the compositions of the fluids and the 
properties of the constituents (e.g., molecular sizes). While NMR data can be used to 
provide detailed information on individual constituents and their properties, most prior art 
NMR data analysis methods only focus on macroscopic properties of earth formation, 
such as where the hydrocarbon and brine zones are, porosity of the earth formations, and 
fractional volumes of pore space; few have focused on a more detailed analysis of 
properties of individual constituents (e.g., molecular size distributions) within a particular 
fluid. 

[0022] As discussed earlier, NMR relaxation rates are dependent on dynamic properties 

of the molecules and the fluids. Thus, NMR data may be used to derive the diffusion 
rates and tumbling rates of the molecules. Because the molecular diffusion rates and 
tumbling rates are sensitive to molecular sizes as well as viscosity of the fluids, NMR 
data may be used to derive information concerning the composition of crude oils in terms 
of molecular sizes. Determining molecular property information from NMR data 
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requires the relaxation time and/or diffusion rate distributions of the hydrocarbon fraction 
of the reservoir and/or borehole fluids. A suitable technique for obtaining these 
distributions in mixtures containing both hydrocarbons and water is the Magnetic 
Resonance Fluid (MRF) characterization method as disclosed in U.S. Patent No. 
6,229,308 Bl issued to Freedman. This patent is assigned to the same assignee hereof 
and is hereby incorporated by reference. 

[0023] The MRF method invokes a Constituent Viscosity Model (CVM), which relates 

relaxation time and diffusion rates to constituent viscosities whose geometric mean is 
identical to the macroscopic fluid viscosity. The validity of the CVM was established by 
Freedman et al using laboratory data acquired on live and dead hydrocarbon mixtures 
and crude oils. These results were reported by Freedman et al in paper number 63214 
entitled "A New NMR Method of Fluid Characterization in Reservoir Rocks: 
Experimental Confirmation and Simulation Results" presented at the 2000 Society of 
Petroleum Engineers Annual Technical Conference and Exhibition meeting. 

[0024] While the MRF method would be used as an example in the following discussion 

of how to derive NMR parameters for individual constituents, one skilled in the art will 
appreciate that other similar methods may be used without departing from the scope of 
the invention. In situations where uncontaminated hydrocarbon samples are available, 
NMR parameters such as relaxation time and diffusion rate distributions may be 
estimated without applying the MRF method. 

[0025] Using the CVM, the MRF method is capable of deriving distribution of 

constituent viscosities in fluid mixtures containing crude oils. The constituent viscosities 
are directly related to the distribution of NMR relaxation times measured on bulk crude 
oil samples and they reflect the complex composition of crude oils as a mixture of many 
different types of hydrocarbon molecules. The use of constituent viscosities simplifies 
the inversion by providing a single set of parameters for characterizing the crude oil 
distributions of bulk relaxation times and diffusion constants. The MRF technique 
provides improvements in all aspects of the standard NMR analysis, including estimates 



7 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2786 



of total porosity, free-fluid and bound-fluid porosity, T2 distributions, and permeability in 
shaly sands. 

[0026] In deriving constituent viscosities, the MRF method uses a general spin-echo 

relaxation model for a formation containing brine, hydrocarbons, and oil-based mud 
filtrate (OBMF). OBMF is commonly used in drilling the borehole. As shown in FIG. 4, 
amplitudes of the spin echoes in an NMR measurement decay exponentially as a function 
of time. Let A? represent the amplitude of the j-th echo acquired during measurement /?. 
Consider the following general relaxation model, 



a j =Z a / ex p 



( j*TE p Y 



l-exp 



_ o 

Yj b k ex p 



k=\ 



j*TE p 



V 



l-exp 
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(1) 



*OBMF 
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f j* TE V 
_ j p 

^ ^2, OBMF if) 



l-exp 



T 

\ 2 \,OBMF J 



where the first, second, and third terms are brine, hydrocarbons, and OBMF signals, 
respectively. This three-phase model explicitly accounts for contributions from 
individual constituents in brine and hydrocarbon phases, but only assumes an average 
relaxation time distribution in the OBMF phase. The simplified term for OBMF is 
justified because experimental measurements in OBMF have shown that NMR relaxation 
time distributions for individual constituents in OBMF are very narrow and can be 
described by a single exponential. 

[0027] The apparent transverse (dipole-dipole) relaxations in any of the three phases 

modeled in Equation (1) include inherent transverse relaxation and relaxation due to self 
diffusion of molecules in the static magnetic field gradient G p . For unrestricted diffusion 
in a uniform magnetic field gradient in the brine, the apparent transverse relaxation rates 
can be written as, 
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i + {yh*g p *te p ) 2 d ^ t) (2) 



KXp) t 2J 12 

Here, T 2 j in the first term on the right-hand side are a set of relaxation times that 
represent the sum of surface and bulk relaxation of the brine phase. The second term is 
the contribution to the relaxation rate from diffusion, where Yh =271*4258 G" 1 s" 1 is the 
proton gyromagnetic ratio and D W (T) is the temperature dependent self diffusion 
coefficient of water in units of cm 2 /s. Note that Equation (2) assumes an unrestricted 
diffusion and a uniform magnetic field gradient G p . One skilled in the art will appreciate 
that corrections to D W (T) for the effects of restricted diffusion and to G p for the effects of 
internal rock gradients can be applied if appropriate. 

[0028] Similarly, the apparent transverse relaxation rates in the native oil (T 2o {7] k ,p)) 

can be written in the form, 

1 >* G :' TE ' y dm p) 



where T 2t0 (tlk) is the bulk relaxation time associated with amplitude b k in the hydrocarbon 
relaxation time distribution, and D o (r\0 is a viscosity dependent diffusion constant. The 
hydrocarbon (crude oil) is usually a non-wetting phase and is not affected by surface 
relaxation. Crude oils are mixtures consisting of many different types of hydrocarbon 
molecules of varying sizes, shapes and molecular weights. See , for example, McCain, 
W. D., The Properties Of Petroleum Fluids, Penn Well Publishing Co., Second Edition, 
Chapter 1, 1990. A molecular-level distribution of constituent viscosities (rj k ) is assumed 
to exist in crude oils. This assumption is based on experimental data that there exists a 
distribution of relaxation times in crude oils. 

[0029] The measured viscosity (r/ 0 ) reflects a macroscopic transport property of the 

crude oil that determines its flow properties and is the quantity that is used in 
hydrodynamic transport equations like the Navier-Stokes equation. Morriss et al showed 
that, for a suite of dead (i.e., not containing dissolved solution gas) crude oils, there 
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exists a strong correlation between the logarithmic mean relaxation times of their 
constituents and the measured viscosities. See Morriss et al. 5 Hydrocarbon Saturation 
And Viscosity Estimation From NMR Logging In The Belridge Diatomite, Paper C 
presented at the 35th Annual Meeting Of The Society Of Professional Well Logging 
Analysis, 1994. The macroscopic viscosity (rj 0 ) of live crude oils is empirically related to 

the logarithmic mean (( T 2, 0 )iogm) of the transverse relaxation time distributions by a 
constitutive equation of the form, 

\ 7 2,Ji ogm 1 



f§\ where a is an empirically determined constituent constant that has been determined by 

~ Morriss, et al to be around 250 (i.e., a = 250 Ks^cp" 1 ), for ( T 2,o)io g m in seconds and rj Q in 

PB 



a f ( ^)R. ") 

centipoise and T the temperature in degrees Kelvin. Thus, c= ^ . The 



empirically derived function f(GOR) accounts for live oils (those containing dissolved 

m 

solution gas) and has been discussed by Freedman et al in paper number 63214 entitled 
p "A New NMR Method of Fluid Characterization in Reservoir Rocks: Experimental 

Confirmation and Simulation Results" presented at the 2000 Society of Petroleum 
Engineers Annual Technical Conference and Exhibition meeting. 

[0030] The //* in Equation (3) are microscopic viscosities that reflect the complex 

composition of crude oils. Analogously with the above equation, the constituent 
viscosities are assumed to be related to the components in the relaxation time distribution 
via the same equation, 

_> = ar h AGO E > 



[0031] The logarithmic mean of relaxation time is defined as, 
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f 2 ,.L„=.0^ S '^-"-» l (6) 
with bk = „ k , where bk are the AT 0 amplitudes in the crude oil bulk relaxation time 

I*. 

Ar=l 

distribution. Substituting Equations (4) and (5) into Equation (6) would yield the 
macroscopic viscosity of the crude oil, rj 0 , which is the logarithmic mean of the individual 
microscopic viscosities, r\k- 

where bk is the "concentration" of the mixture constituent with viscosity The 
macroscopic viscosity rj 0 is similar to the high temperature limit for the viscosity of a 
mixture according to the "Arrenhius mixing rule" see A. Bondi, Physical Properties of 
Molecular Crystals, Liquids, and Glasses, pp. 348-349, 1968. 

[0032] The dependence of the relaxation times on viscosity and temperature in Equations 

(4) and (5) is consistent with the experimental observations and theoretical predictions of 
Bloembergen, Purcell, and Pound, Relaxation Effects in Nuclear Magnetic Resonance 
Absorption, Physical Review, vol. 73, no. 7, pp. 679-712, 1948. 

[0033] Stokes-Einstein diffusion theory predicts that diffusivity is related to temperature 

and viscosity according to the equation: D = kT/6nr\R, where k is the Boltzmann 
constant, R is the radius of the spherical particle, and T is the temperature in degrees 
Kelvin. Similar to the Stokes-Einstein equation, the self-diffusion constants for the crude 
oils, D 0y and for constituents in the crude oil, D 0 (rj k ), are assumed to have the same 
dependence on Tlrj^ Therefore, for crude oils, 

D = ^xlO" 5 (7) 

Vo 

where b is a constant, D Q is the measured crude oil diffusion constant in cm /s, and T is 
the temperature in degrees Kelvin. The empirical constitutive constant for crude oils , 
b = 5.05 x\0~ 3 cm 2 s~ l cpK~ l , is given by Freedman et al in paper number 63214 entitled 
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"A New NMR Method of Fluid Characterization in Reservoir Rocks: Experimental 
Confirmation and Simulation Results/' presented at the 2000 Society of Petroleum 
Engineers Annual Technical Conference and Exhibition meeting. 

[0034] Analogously to the relationship between the macroscopic diffusion constant (D 0 ) 

and the macroscopic viscosity (rj 0 ), the microscopic constituent diffusion constants of 
crude oil mixtures are related to the microscopic constituent viscosities (effective 
viscosities) according to the following equation: 

DM = — xlO" 5 (8) 

Equation (8) implies that there exists at the molecular level a distribution of diffusion 
constants in the crude oil mixture analogous to the distribution of relaxation times. 
These distributions of diffusion constants and relaxation times can be derived from the 
measured NMR data by iteratively fitting a model as shown in Equation (1) to these data 
using the method disclosed in U.S. Patent No. 6,229,308 Bl. 

[0035] The above-described MRF method is just one way to obtain the distributions of 

the diffusion constants and the relaxation times. The MRF method is particularly 
appropriate when NMR data are obtained from mixed fluids (e.g., water, drilling fluid 
filtrates, and oil). If the oil sample is not contaminated with other fluids (e.g., crude oils), 
there would be no need to apply the MRF method to obtain the distributions of diffusion 
constants and the relaxation times. 

[0036] Once the distribution of the diffusion constants and the relaxation times are 

estimated, they can be used to further derive the molecular properties of the individual 
constituents according to embodiments of the present invention. Molecular properties as 
used herein refer to molecular size, carbon number, and weight, i.e., those properties 
related to physical dimensions of the molecules. For example, assuming a spherical 
molecule with a radius of R, the individual diffusion constant is related to the radius R 
(hence, the molecular size) according to the Stokes-Einstein equation: D 0 = kT/67trjR, 
where D 0 is the observed diffusion constant, k is the Boltzmann constant, T is the 
temperature in degrees Kelvin, and rj is the viscosity in centipoise. 
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[0037] Therefore, to derive molecular size information from the NMR data, the 

relaxation times (T2 k ) and diffusion rates (D k ) of individual molecules (indexed by k) in a 
hydrocarbon mixture can be approximated as 

T2 = e^Il (9) 

n a f(GOR)N? 



D k = J^l (10) 



where Nk is the number of carbon atoms in the kth constituent, a '(T) and b '(T) are 
functions of temperature, T, rj is the fluid viscosity, and of, ft and ff are (as yet) 

□ unknown exponents. Note that these expressions are generalizations of the Stokes- 

O 

ry Einstein and Bloembergen relations for diffusion and spin relaxation of spherical 

5 particles in liquids. Equation (9) includes a factor that represents the gas/oil ratio (GOR). 

yg This is included because it is known that GOR is an important parameter in determining 

S3 

7 the relaxation time dependence on viscosity and temperature. See Lo et al., Relaxation 

^ Time And Diffusion Measurements of Methane And N-Decane Mixtures, The Log 

ru 

Analyst, pp. 43-46, November-December, 1998; see also U.S. Patent No. 6,229,308 Bl. 
% Within the context of the CVM approach, the denominators of Equations (9) and (10) are 

H 5 proportional to constituent viscosities. 

[0038] With the approximations in Equation (9) and (10), it is tempting to assume that 

the exponents a 9 and <p' are equal to 1 and then make some correlation between Nk and 
molecular "radius," R, to mimic the ideal spherical particle expressions. Such 
assumption suggests a dependence on the inverse of rjR for D and on the inverse of r/R 3 
for T2. However, this approach immediately fails for mixtures because it implies that the 
geometric mean of the Nk distribution is independent of the details of the mixture. This 
contradiction follows directly when the expressions in Equations (4) and (7) are 
respectively substituted into the constituent Equations (9) and (10). 

[0039] By incorporating empirical results concerning the relation between viscosity and 

the geometric mean relaxation times and diffusion rates, i.e., T2 LM and D LM are linearly 
dependent on T/>/, Equations (9) and (10) can be reformulated as, 
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= a(T)T2 a LM [f(GOR)r (n) 

N, = (12) 

Note that Equations (11) and (12) are obtained from Equations (9) and (10) simply by 
introducing the empirical expressions, T2lm = a 'T/rj f(GOR) and Dlm = b T/7, re-ordering 
the variables, and defining a new set of exponents (unprimed). 

[0040] Equations (11) and (12) indicate that if the exponents, a, /?, fa and 0, are known, 

then Nk (carbon number) may be obtained from the geometric means of T2 and D (i.e., 

T2 L m and D L m)- Because relationships shown in Equations (11) and (12) are not 
dependent on the exact natures of the individual constituents in the mixture, it should be 
possible to derive these exponents, , a; /?, fa and 0, using a simple model mixture system. 
Once these exponents are derived, they may be used in other systems of similar 
compositions. 

[0041] As an example, rough estimates for the exponents, a, /?, fa and 0, which would be 

useful in other hydrocarbon mixture systems, can be derived from an experimental 
mixture of squalene (C30) and hexane (C6). For this mixture, GOR = 0 and f(GOR) = 1. 
Data for the pure squalene (C30) and hexane (C6) and three mixtures of these two 
components at 30° C are given in Table 1. Note that any suite of hydrocarbon samples, 
including crude oils, may be used to derive such parameters. 

Table 1. 



Relaxation Time and Diffusion Rates for the Hexane-Squalene System 



Fluid 

(% C6) 


T2 C 30 
(ms) 


T2 C6 
(ms) 


DC30 

(cmV) 


D C6 
(cmV 1 ) 


T2lm 
(ms) 


Dlm 
(cmV) 


n(cp) 


0 


300 




1.0x10"" 




300 


i.oxio- b 


11 


38 


900 


4600 


3.7xl0" 6 


1.3xl0" 5 


1670 


6.0x1 0 6 


1.7 


50 


1180 


5800 


4.9x1 0" 6 


1.7xl0' 5 


2610 


9.1xl0 -6 


1.2 


69 


1690 


7500 


9.9xl0 _6 


3.0x1 0' 5 


4780 


2.7X10" 5 


0.6 
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100 



9800 - 4.6x1 0' 5 9800 " 4.6x1 0" 3 0.3 



[0042] Once, the T2, T2 LM , D, D L m values for squalene (C30) and hexane (C6) are 

available, they can be plugged into Equations (11) and (12) to produce the rough 
estimates of the exponents, , a, ft, <f>, and 0. From the dataset shown in Table 1, suitable 
exponents and pre-multiplying functions are determined, and the expressions for N k (at 
30° C) as shown in Equations (1 1) and (12) are simplified as: 

N t - 650 T2 ™ (13) 

T2 k . 



§ N m (H) 

% [0043] Equations (13) and (14) can be used to estimate molecular sizes of individual 

00 constituents in a similar mixture. Using Equations (13) and (14) and the measured 

relaxation and diffusion data shown in Table 1, the N k values for the squalene-hexane 
system are given in Table 2. 

Table 2 



Number of Carbons Estimated for Squalene (C30) and Hexane (C6) from Relaxation 
Time (T2) and Diffusion Rate (D) Measurements 



Fluid (%C6) 


N C30 (T2) 


N C 3o(D) 


N C6 (T2) 


N C6 (D) 


0 


37 


40 






38 


30 


34 


5.8 


5.1 


50 


28 


34 


5.7 


5.2 


69 


27 


34 


6.0 


6.6 


100 






6.6 


5.9 



[0044] As shown in Table 2, reasonable values for Nk can be obtained using either 

relaxation times or diffusion rates. It should be noted that Equations (13) and (14) and 
the example used here are merely intended to illustrate the basic concept of how to obtain 
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these values; they should not limit the invention. One skilled in the art will appreciate 
that the derived exponents in the above example were rounded to give simple half-integer 
fractions for simplicity. Alternative optimized values may be derived which should 
provide better estimates for N*. However, as shown in the above example, even with the 
simplified rough estimates, the Nk values can be obtained with reasonable accuracy. 
Therefore, the rough estimate approach as shown here should be sufficient for most 
situations. 

[0045] The pre-multipliers, 650 and 0.04, in Equations (13) and (14), respectively, in 

principle are valid only at the measurement temperature of 30° C. However, the 
temperature dependence of these values is relatively weak (approximately T with T in 
degrees Kelvin). Thus, these values may be used in a temperature range around 30° C. 
For crude oils, the optimum exponents and pre-multipliers might differ slightly from 
those derived for simple bi-component mixtures. Equations (13) and (14) are not strictly 
compatible with the CVM equations because T2 k and D k depend differently on Nk and r\. 
The diffusion constants Dk 9 as shown in Equation (14), exhibit a weaker dependency on 
Nk than do the relaxation times, as shown in Equation (13). This result agrees 
qualitatively with the ideal spherical particle relations. In view of the many 
approximations and assumptions implicit in this kind of model, any resulting "carbon 
number" (or molecular size) distribution should probably be regarded as an approximate 
indicator rather than a definitive and accurate breakdown of molecular composition. 

[0046] Equations (9) - (12) should be regarded as particular implementations of the 

method. Alternative mathematical expression relating relaxation times and diffusion 
rates with molecular size, carbon number or other constituent property could also be 
derived and calibrated to hydrocarbon mixtures. It is also feasible to determine molecular 
properties from measured NMR data using model independent pattern recognition 
methods such as neural networks. In this approach, there are no model dependent 
equations (e.g., equations 13 and 14). Instead a "training data set" of molecular 
properties versus NMR and diffusion properties is used to train a neural network to 
predict molecular properties given NMR data on a sample outside of the training set. 
Any commercially available neural network software (such as that available from the 
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Mathworks, Inc. at www.mathworks.com) may be adapted to determine molecular 
properties without the need to invoke model equations. These methods might also 
incorporate additional data derived from other measurements, for example NMR 
spectroscopy or optical analysis. 

[0047] FIG. 5 summarizes steps involved in methods for evaluating molecular size 

distribution in a hydrocarbon mixture according to embodiments of the invention. First, 
NMR data are collected (process 41 in FIG. 5). This may be performed with any suitable 
NMR logging tool. For example, embodiments of the invention are applicable to 
wireline and while-drilling NMR tools as well as any suitable NMR module sampling 
tool such as the Modular Formation Dynamics Tester (MDT™ from Schlumberger 
Technology, Houston, TX). In fact, the absence of fluid-rock interactions, which 
complicate the MRF analysis for wireline or LWD-NMR tools, should lead to greater 
accuracy and robustness of the method in a sampling module or laboratory analysis 

[0048] Once the NMR data are collected, they are analyzed using an inversion method to 

derive individual constituent dynamic parameters (e.g., Tl, T2, T1/T2, and diffusion 
constants; process 42 in FIG. 5). As discussed earlier, the MRF method or any similar 
method may be used for this purpose. Note that the MRF technique and the extension to 
it described herein are able to provide real-time information on reservoir fluids (e.g., 
viscosity, molecular composition) that at present can only be provided by lengthy 
pressure-volume-temperature (PVT) analysis performed in laboratories. 

[0049] Finally, the individual constituent dynamic parameters (e.g., Tl, T2, T1/T2, and 

diffusion constants) may be used to derive the molecular size information (process 43 in 
FIG. 5). As discussed above, the molecular sizes can be correlated with the transverse 
relaxation times and the diffusion constants according to Equations (11) and (12). The 
exponents in these equations can be estimated using a model mixture having similar 
components and/or properties {i.e., hydrocarbons) under similar conditions (e.g., 
temperature). Rough estimates of these components ("empirical parameters") would be 
sufficient. Having these exponents, Equations (11) and (12) may be simplified to those 
like Equations (13) and (14). The individual constituent dynamic parameters (e.g., Tl, 
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T2, T1/T2, and diffusion constants) in the mixture of interest derived from process 42 
may then be used to calculate the molecular sizes of constituents (or their distribution in 
the mixture). 

[0050] FIGs. 6A - 6D illustrate results obtained with methods of the present invention as 

compared with those obtained with gas phase chromatography (GPC) that is commonly 
used in a PVT laboratory for such analysis. This example also illustrates that, in addition 
to wellsite applications, methods of the invention are equally applicable to standard 
laboratory analysis of crude oil samples. 

[0051] FIG. 6 A shows the T2 relaxation time distributions of three different dead crude 

oil samples. Curves Al - A3, representing samples 1 - 3, have geometric mean T2 L m 
values of 105 ms, 239 ms, and 603 ms, respectively. These measurements were obtained 
at 300° K. 

[0052] Curves Al - A3 were analyzed with methods of the invention and the resultant 

molecular weight distributions are shown as curves B2 - D2 in FIGs. 6B - 6D, 
respectively. These results were obtained using a(T) = 170, T= 300° K, a= 0.51, and /? = 
0.90 in Equation (11). For comparison, the corresponding molecular weight distributions 
as obtained with GPC analysis are superimposed in FIGs. 6B - 6D as curves Bl - Dl, 
respectively. It is clear from these figures that the NMR methods according to 
embodiments of the present invention produce results similar to those obtained from the 
GPC method. It should be noted that a GPC analysis takes several hours, while an NMR 
measurement takes only a few minutes. Furthermore, GPC analysis relies on correlation 
of component retention times with their molecular weight, which may be no more 
reliable than methods of the present invention. The NMR approach also has further 
advantages that it is non-destructive, and sample handling is relatively simple. 

[0053] Although the example in FIG. 6 uses NMR transverse relaxation time (T2), one 

skilled in the art would appreciate that other NMR parameters (Tl, T1/T2, or diffusion 
constant) may be used. Furthermore, while this example used a mixture of hydrocarbons 
(crude oils), methods of the invention may also be applied to other liquid mixtures. 
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[0054] While the invention has been described using limited examples, those skilled in 

the art, having the benefit of this disclosure, will appreciate that other methods can be 
devised without departing from the scope of the invention as disclosed herein. 
Accordingly, the scope of the invention should be limited only by the attached claims. 
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